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Abstract 

RIGl and MDA5 have emerged as important intracellular innate pattern recognition receptors that recognize viral RNA and 
mediate cellular signals controlling Type I interferon (IFN-I) response. Buffalo RIGl and IV1DA5 genes were investigated to 
understand the mechanism of receptor induced antiviral response. Sequence analysis revealed that RIGl and MDA5 
maintain a domain arrangement that is common in mammals. Critical binding site residues of the receptors are evolutionary 
conserved among mammals. IVlolecular dynamics simulations suggested that RIGl and IVIDAS follow a similar, if not 
identical, dsRNA binding pattern that has been previously reported in human. IVloreover, binding free energy calculation 
revealed that IVIDAS had a greater affinity towards dsRNA compared to RIGl. Constitutive expressions of RLR genes were 
ubiquitous in different tissues without being specific to immune organs. Poly l:C stimulation induced elevated expressions 
of IFN-p and IFN-stimulated genes (ISGs) through interferon regulatory factors (IRFs) mediated pathway in buffalo foetal 
fibroblast cells. The present study provides crucial insights into the structure and function of RIGl and IVIDAS receptors in 
buffalo. 
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Introduction 

Establishment of an immune response begins with the 
discrimination between "seLP' and "non self components within 
an organism. To recognize "non seU" pathogen-associated 
molecular patterns (PAMPs), vertebrates have evolved a number 
of germ-line-encoded receptors called pattern recognition recep- 
tors (PRRs). Over the past decades, several types of PRRs have 
been identified, the recent being cytosohc retinoic acid inducible 
gene 1 (RIGl)-like helicases (RLHs) including RIGl and melano- 
ma differentiation-associated gene 5 {MDA5) [1]. RIGl and MDA5 
along with another gene IGP2 (Laboratory of Genetics and 
Physiology 2) constitutes cytosolic RIGl-lilie receptors (RLR) 
family, capable of sensing double stranded RNA (dsRNA) 
signatures of viruses. RIGl and MDA5 share sequence homology 
and transduce downstream signal through a common adaptor 
called mitochondrial antiviral signaling protein (MAVS) that 
activates interferon (INF) signaling pathways [2,3,4]. RLRs differ 
from other classes of PRRs viz- Toll-like receptors (TLRs) and 
NOD-hke receptors (NLRs), in that these receptors lack repetitive 
receptor elements such as leucine rich repeats (LRRs) as 
recognition platform of PAMPs. Instead, they recognize viral 
dsRNA with a central DExD/H box helicase domain which is 
characterized by seven conserved motifs (I/Walker A, la, II/ 
Walker B, III IV, V and VI) [5,6]. This conserved "helicase" core 



of RIGl and MDA5 is flanked by two tandemly arranged N- 
terminal caspase activation and recruitment domains (CARDs), 
and a cysteine-rich C-terminal regulatory domain (CTD) [7]. 
LGP2 has similar domain architecture as RIGl and MDA5, but it 
lacks N-terminal CARDs [8] . In normal cell, CTD plays a role in 
autoinhibition of RIGl and MDA5 through intramolecular 
interactions that mask the CARD signaling [9,10]. Study on 
crystal structures of CTD suggests that it plays a prominent role in 
high affmity binding and selectivity for dsRNA [8,11]. However, 
MDA5's CTD displays lower magnitude of affinity towards 
dsRNA than that of fuU-length MDA5 [12]. Thus, despite 
homology, they differ in molecular mechanisms for dsRNA 
recognition and play non-redundant function in antiviral immu- 
nity [13,14]. MDA5 detects long duplex RNAs in the genome of 
double-stranded RNA (dsRNA) viruses or dsRNA replication 
intermediates of positive-strand viruses, e.g. Encephalomyocarditis 
and vesicular stomatitis viruses [15,16]. In contrast, RIGl detects 
the 5' triphosphate group (5'ppp) and blunt end of short dsRNA 
or single strand RNA (ssRNA) hairpins, which are often present in 
a wide range of negative strand viruses, {e.g. influenza A virus) and 
in some positive and double strand RNA viruses [17,18]. The role 
of LGP2 in antiviral response is less clear as compared to RIG 1 
and MDA5. Since LGP2 lacks N-terminal CARDs, it cannot 
mediate MAVS based signal transduction. Nonetheless, it is 
considered to act as a suppressor of RIGl and MDA5 elicited 
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signaling [19,20]. Other studies, however, have shown LGP2 as a 
positive inducer of RIG I and MDA5 amplifying their responses 
via a positive feedback mechanism [21,22]. 

Buffalo (Bubalus hubalis) Ls an important milk and meat 
producing animal in many Mediterranean and East Asian 
countries including India. Besides, Indian buffaloes are well- 
known for their exceptional ability of disease resistance [23,24]. 
However, the genetic basis of disease resistance is still elusive due 
to paucity of information on molecular mechanism(s) of host- 
pathogen interaction. Here, we studied the mechanism of viral 
RNA recognition by buffalo RIGl and MDA5 receptors through 
molecular dynamics (MD) simulations. We investigated the 
expression profile of c:orresponding genes constitutively, and in 
response to a stimulus in order to understand the mechanism of 
RLR mediated antiviral response in buffalo. 

Materials and Methods 

Materials 

Taq DNA polymerase, lOX buffer, dNTP were purchased from 
New England Biolab (MA, USA) and Fermentas Maxima SYBR 
Green qPCR Master Mix (2X) was obtained from Thermo Fisher 
Scientific Inc. (PA, USA). Unlabelled goat anti-bovine polyclonal 
antibodies against RIGl and MDA5 were obtained from Santa 
Cruz Biotechnology Inc, (CA, USA) and horseradish peroxidase- 
labeled rabbit anti-goat antibody was obtained from Sigma- 
Aldrich (St. Louis, MO, USA). Poly (I:C)/LyoVec HMW 
(endotoxin leveKO.OOl EU/^g) was obtained from Invivogen 
(San Diego, CA). Dulbecco Modified Eagle's medium (DMEM), 
Trypsin-EDTA, foetal bovine serum (FBS) and Epidermal Growth 
Factor (EGF) were obtained from Sigma-Aldrich (St. Louis, MO, 
USA). L-Glutamine (Glutamax lOOx) was purchased from 
Invitrogen corp., (Carlsbad, CA, USA). PeniciUin-G and strepto- 
mycin were obtained from Amresco (Solon, OH, USA). Filters 
(Millex GV. 0.22 |xm) were purchased from Millipore Pvt. Ltd., 
(BUlerica, MA, USA). AU other reagents were of analytical grade. 

PCR Amplification 

An PGR amplifications were performed in 25 |a,l reaction 
volume. Each reaction contained 2.5 |il lOX buffer, 200 |.lM of 
dNTPs, 0.5 |J,1 of each primers (10 pmol), 0.5 units of Taq DNA 
polymerase and nuclease free water to bring the total volume to 
25 Around 100 ng of cDNA was used as template. The details 
of primers for different fragments of RIGl and MDA5 genes are 
summarized in Table 1. Thermal cycling was performed on an 
ABI Veriti PCR machine where a touchdown protocol with 
annealing 58°C was followed (Table SI in File SI). The PCR 
products were resolved on a 1.5% agarose gel (Figure SI in File 
SI). The PCR products were cloned on to pTZ57R/T vector, 
plasmids were screened using PCR, and positive plasmids were 
custom sequenced. 

Immunohistochemistry (IHC) 

Spleen, tonsil and thyroid samples were collected from 
Municipal Corporation of Delhi (MCD) Slaughter House, 
Ghazipur, New Delhi, India with permission for research use. 

The samples were fixed with 10"/o ncutral-buffc'rc'd formalin and 
embedded in paraffin. For antigen retrieval, thin sections of tissues 
(~3 |J.m) wi;r(' deparaffinized with xylene and rehydrated by 
washing with ethanol (95%) and MiliQ_ water. Non-specific sites 
were blocked by overnight incubation of sections with a blocking 
reagent (1% BSA, 10% FBS-PBS) at 4°C. Tissue sections were 
then stained with either anti-RIGl (1:100 dilution) or anti-MDA5 
(1:150 dilution) for 1 h at room temperature. Slides were washed 



three times with Tris-buffered saline (pH 7.6) supplemented with 
0.05% Tween-20. Thereafter, the shdes were treated with 
horseradish peroxidase-labeled rabbit anti-goat antibody for 1 h 

and washed as described above. The slides were then mounted 
with DAB (3, 3'-diaminobenzidine) substrate chromogen solution 
for 5 minutes. Counter-staining was performed with Mayer's 
haematoxylin. The slides were rinsed with water, dried and 
visualized. 

Poly l:C Treatment of Foetal Fibroblasts 

Gravid uteri (2-3 months of pregnancy) were obtained from 
abattoir house. New Delhi Municipal Corporation. For poly I:C 
treatment, buffalo foetal fibroblast cells at about 70% confluence 
were harvested and seeded at a concentration of 1 xlO^ceUs/ml 
into four well culture plates (Text SI in File SI). The cells were 
cultured in DMEM (10% FBS) for 24 h at 37°C. After 24 h, 
serum containing medium was decanted and cells were washed 
twice with Dulbecco's phosphate buffer saline (DPBS) solution. 
Cells were then maintained with serum free DMEM medium at 
37°C for 6 h. The medium was subsequently discarded and cells 
were maintained either with serum free DMEM medium (control 
group) or serum free DMEM medium containing 2 |ig/ ml of Poly 
(I:C)/LyoVec (treatment group). The experiment was performed 
with four replicates for each treatment or control group. 

Real Time PCR (qRT-PCR) 

Real time PCR primers were designed by aligning gene 
sequences of several mammals including cow, pig, mouse, buffalo, 
and human (Table 2). Equal amount of RNA (quantified by Qubit 
fluorometer, Invitrogen), wherever applicable, were used for 
cDNA preparation (Superscript III cDNA synthesis kit; Invitro- 
gen). All qRT-PCR reactions were conducted on a Light Cycler 
480 II Real-Time PCR machine (Roche Diagnostics, USA). Each 
reaction consisted of 2 |J,1 cDNA template, 5 \A of 2X SYBR 
Green PCR Master Mix, 0.25 |il cadi of forward and reverse 
primers (10 pmol/ H-l) and nuclease free water for a final volume of 
10 nl. Each sample was run in duplicate. Analysis of real-time 
PCR (qRT-PCR) was performed by delta-delta-Ct (AACt) method 
[25]. 

Molecular Dynamics Simulations 

The initial structures of RIGl and MDA5 were obtained using 

homology modeling (Text S2 in File SI). Molecular dynamics 
(MD) simulations of modeled RIG 1 and MD A5 were carried out 
in GROMACS 4.5.5 [26] using Charmm27 all-atom force field. 
All the structures were hydrated using SPC216 water molecules in 
two separate cubic boxes with a distance of 10 A between the 
proteins and box edges. Physiological ionic strengths (0.15 M) of 
counter ions were added to neutralize total charge of the 
simulation systems. The electroneutralized systems were subjected 
to steepest descent energy minimization to remove steric conflicts 
between atoms. The minimized models were then subjected to 
MD simulations using NPT ensemble with constraints on the 
backbone heavy atoms for 100 ps, keeping the number of particles 
(N), the system pressure (P), and the temperature (T) constant. The 
temperature of the systems was coupled to 300 K via Berendsen 
temperature coupling scheme, and the pressure of the systems 
were maintained at 1 bar using ParrineUo-Rahman algorithm. 
Finally, the models of RIGl, and MDA5 were subjected to 2 ns 
MD simulations without any constraints on the backbone atoms. 
Trajectory and subsequent data analyses were performed using 
VMD 1.9.1 [27], Grace 5.1.23 (http://plasma-gate.weizmann.ac. 
il/ Grace/), and Pymoll.3rc2 (The PyMOL Molecular Graphics 
System, Version 1.3 Schrodinger, LLC) software. Principal 
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Table 1. Primers used for amplification of buffalo RIG1 and MDA5 genes. 





Gene/Fragment 


Primer 


Sequences (5 -3 ) 


Length (bp) 


T„ (°C) 


RIGIFragmentI 


Forward 


CAGGCATGACGGCAGAGCAGCG 


22 


62 




Reverse 


TTGATAATGAGGGCATCATTGTATTTCCGCA 


31 


59 


RIG1Fragment2 


Forward 


GTGTTTCAGATGCCAGACAAAGAGGAAGA 


29 


60 




Reverse 


TCTTCCTCTGCCTCTGGTCTGGATCAT 


27 


61 


RIGl Fragments 


Forward 


GTCGCCGATGAAGGCATTGACATTGC 


26 


61 




Reverse 


CCTGAGCCCAAGGGGACATTTCTGC 


25 


63 


MDA5 Fragmentl 


Forward 


AGATGTCGTCGGACGGGAGTTCAACAGA 


28 


63 




Reverse 


CTGAATCACTTCCCATGGTGCCTGAATC 


28 


61 


MDA5 Fragment2 


Forward 


GTTTGGCAGAAGGAAGTGTCAGCTGCT 


27 


61 




Reverse 


CCTCATTGTACTTCCTCAGATGTTCTGCAC 


30 


62 


MDA5 Fragments 


Forward 


AGCTGCAAAAGAAGGAAATCGCAAAGATCGTG 


33 


62 




Reverse 


GAAATCTTGAATCAAATGTCAGTCCTCATCAC 


32 


59 



doi:l 0.1 371 /journal.pone.0089788.t001 



component analysis of MD trajectories was performed using 
'g_covar' and 'g_anaeig' programs of GROMACS (Text S2 in File 
SI). 

Binding Free Energy Calculation 

The binding free energies (AGi,ind) of RIGl-dsRNA and 
MDA5-dsRNA were calculated with GMXAPBS tool [28]. The 
detailed procedure is outlined in Text S2 in File SI. 



Statistical Analysis 

All statistical analyses were carried out using SYSTAT 13.1 
software (SYSTAT Software Inc.) Analysis of variance (ANOVA) 
was used to test between groups and hour intervals. Fischer's 
restricted least significant differences criterion was used to 
maintain the a priori type I error rate of 0.05. 



Table 2. Primers used for mRNA quantitation of different genes by Real time PGR (qRT-PCR). 





Gene 


Primer 


Sequence (5 -3 ) 


T„, CC) 


Product Size (bp) 


RIGl 


Forward 


CTTGCAAGAGGAATACCACTTAAACCCAGAGAC 


63 


150 




Reverse 


TTCTGCCACGTCCAGTCAATATGCCAGGTTT 


63 




MDA5 


Forward 


TCTGCTTATCGCTACCACAGTGGCAGA 


61 


135 




Reverse 


TGCTCTCATCAGCTCTGGCTCGACC 


63 




LGP2 


Forward 


CCACCACGTCAATGTGAACCCCAACTTC 


63 


176 




Reverse 


TGAGCACTGGCAGCTTCACTGACTTGTAGAT 


63 




MAVS 


Forward 


CATCAGGAGCAAGACACAGAACTGGGCA 


63 


115 




Reverse 


AACGGGCCAAGGGCTGGAAGGAGACAC 


66 




IRF3 


Forward 


GATACTGCCCTGGCTGATATCTCAGCTG 


63 


235 




Reverse 


GTTCAGGGCAGACCGGAAATTCCTCTTC 


63 




IRF7 


Forward 


GCCTCCTGGAAAACCAACTTCCGCTG 


63 


340 




Reverse 


GCAGATGGTCCTCCAAGCAGCTCTG 


63 




IFNfi 


Forward 


CACCACAGCTCTTTCCAGGAGCTACA 


61 


168 




Reverse 


GGAACTGCTGTTCTTGCTTCATCTCCTC 


61 




ISG15 


Forward 


AGATCAATGTGCCTGCTTTCCAGCAGCG 


63 


162 




Reverse 


GACCCTTGTCGTTCCTCACCAGGATG 


61 




ISG54 


Forward 


AAGTGCACGGCAATCATGAGTGAGACCA 


61 


94 




Reverse 


CCTCTACCAAGTTCCAGGTGAAATGGC 


61 




ISG56 


Forward 


CTATGTGAAACACCTGAGAGGCCAGAATGA 


62 


132 




Reverse 


CAGGCATAGTTGCCCCAGGTAACCAG 


63 





doi:1 0.1 371 /journal.pone.0a89788.ta02 
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Results 

Sequence Analysis of Buffalo RIG1 and MDA5 Genes 

The amino acid sequences of bufTalo RIG 1 and MDA5 were 
translated from respective nucleotide sequences identified in this 
study (GenBank Accession ID: KF5 17376 and KF5 17377; Text 
S3 in File SI). The putative conserved domains and critical 
binding site residues within RIGl and MDA5 were identified 
using DELTA BLAST and Conserved Domain Database (ODD) 
of NCBI. Comparison of amino acid conservation among difiFerent 
domains of buffalo RIGl and MDA5 sequences revealed that 
CARD, helicase, and CTDs of these genes share 31.1, 34.9 and 
21.5% homology, respectively (Figure 1). Putative conserved 
domains along with conservation of critical binding site residues 
among buffalo RIGl, MDA5 and different model organism has 
been shown in Figure lA. From Figure lA it is clear that critical 
binding site residues of both RIGl and MDA5 are conserved 
across mammals. Moreover, the residues responsible for ATP, 
Mg"*"*", and RNA binding are almost identical. Phylogenetic trees 
constructed by maximum likelihood method (MEGA5.2 software) 
[29] revealed a similar pattern of evolution for RIGl and MDA5 
in different model organisms except perrisodactyla (horse), where, 
RIG 1 was clustered with cetartiodactyla and MDA5 with primates 
(Figure IB). 



Structural Analysis of RIGl and MDA5 Proteins 

The modeled buffalo RIGl and MDA5 exliibits characteristic 
tertiary folds containing HEL-1, HEL-2, HEL-2i, and the CTD 
(Figure S2 in File SI). The stereochemical qualities of the modeled 
proteins were satisfactory, indicating that the models could be used 
in simulation studies (Figure S3 in File SI; Table S3 in File SI). 

MD simulations were performed on free and dsRNA bound 
RIGl and MDA5 to analyze their structural organization and 
intermolecular interactions. Analysis of MD trajectories revealed 
that the backbone as well as the individual residues of RNA bound 
proteins fluctuate less compared to free proteins (Figure S4 in File 
SI). This observation prompted us to study the global motions of 
individual domains during MD simulation. Principal component 
analysis (PC A) revealed that CTD of free RIGl shows maximum 
displacements (Figure 2A). However, CTD of RNA bound RIGl 
was comparatively rigid. Further, it was found that HEL-2i 
domain of both free and RNA bound proteins display movements 
(Figure 2A). As can be seen from Figure 3B, all four domains of 
RIGl, i.e. HEL-1, HEL-2, HEL-2i, and CTD showed rotations 
and translations upon RNA binding. Similarly, free as well as 
RNA bound MDA5 showed noticeable motions of individual 
domain (Figure 2B). In addition, the V shaped linker of RNA 
bound proteins showed major displacements as compared to 
unbound proteins. Collectively, all four domains were found to be 
involved in dsRNA recognition. 



Cnttle 

Mouse 

Human 

Rhesus 

Cat 

Sheep 

Pie 

Zebra Fine 
Xeiiopus 



DEAD BOX HELICASE 



T" 



CARD1 C«R0-2 

'T ' 



11-198 



322 - 828 




.Y SSJIIAVC AC, 

.Y SWIIAVQ AC. 

.Y ES>JIVGVIi AC. 

.Y C^nVTVC AC. 

.Y fiWIVTVC AC. 

-Y aWIVTVC AC. 

.Y ei-rilAVQ AC. 

.Y QWIrWC AC. 

.Y awmvQ EC. 

.Y KVJIVSTQ KC. 

.Y aWIVDTg KC. 




KNKEMRliS SLPTYT 



-I,i1!:edar: ' 
..ilti edar- 
..vlii edaq' 1 

..iLII.-EHARI,- 
.ALII3EHABLF 
-.iillSEHAK-lr- 

..'U.II.yEHRBp'tt. 
.ALIIMEHABl 
.ALllSEDAEll 
. SIMrNDDARTI-i 




Buffalo RIG-1 
933 AA (Partial) 



Buffalo MDA-S 
1021 AA 



Kadei I. 

KADEI I. 

KADEI I. 

«ADEI I. 

KADEI I . 

KAnEI I. 

KADEI I. 

K^DEI I. 

KADEI I. 

K/iDEI I. 

^jADEI I. 



TRG 3 AR 

TRG 2 AR 

TRG 3 AR 

TRG 3 ^R 

TRG 3 AR 

TRG S ?iR 

TRG 3 

TRG 3 AR 

TRG 3 AR 

MI VRG 3 ?iK 

Ml|vRG 3 jy^ 



B 




ATP binding Mg++ RNA binding Regulatory domain 

binding site Interface 



Helic 



e domain interface Nucleotide binding site ATP binding site 




Human 
C3iimp 
Pig 
Cat 


PTsic 
PTa 
Pl9 

pt9 

pt9 


Hrv. 
Hrv. 
Hrv. 
Hrv. 




BT. .Y 


«T. .Y 
HT. -Y 
HT. . Y 


Buffalo 




Cattle 


PIW 


Hrv. 




HT. -Y 
BT. . Y 


Sheep 
Horse 


PT« 
PT« 


Hrv. 
Hrv. 




HT. .Y 
HT. .Y 


Mouse 


PTsI 


Hrv. 


HT. .Y 


Chickeu 
Xenopus 


ptII 

PTia 


Hrv. 
Brv. 




HT. . Y 
HT. -Y 


Conservation 






i 





^Jl'QME KA. .A EIHKE VC. .AL| 
KA. .A bHkE VC 

X,IOME i<A- .A dHkD VC 

rrjiQME KA-,A eHkd VC 

X,IQME fCA. .A eHkD VC 

IJQME E^. .A eHkD VC 
L.IQME .A eHkD VC 

IWAJfME KA- .A eBkD VC 

[vTAigME [-CA. .A dHkD VC 

1REE FA.- .A EiSke|vC 

1QTD TA- .a|ek|kEHVC 




AYT LE. .F' 

AYT LE. 

IRI AYN LK. 

|tI AYN LE. .F' 

AYN LE. 

AYN LE. 

AYN LE. 

AYN LE. .F' 

AYS LE. 

AYN LN. .F' 

SLI LR 




• GL IV JUIG R aD 

•GL IV ARG A 

.GL IV ARG A 

.Gl IV ARG A 

.GL IV ARG A 

-GL IV ARG A 

.GL IV ARG A AD 

.GL IV ARG A AD 

.GL IV ARG A 

.GL IV ARG A AD 
.GL Iv|aRC^ ad 




Mg++ 
binding site 



Nucleotide binding 



Figure 1 . Amino acid sequence analysis of RIG1 and MDA5 receptors. (A) Sequence comparison of functional domains and binding residues 
of RIGl and iVIDAS. Multiple sequence alignment was performed by MAFFT web server and binding site residues (highlighted) were identified by 
DELTA-BLAST and CDD database of NCBI. (B) The evolutionary history was inferred by using the Maximum Likelihood method of MEGA5 based on the 
JTT matrix-based model. The trees with the highest log likelihood are shown. Initial tree(s) for the heuristic search were obtained automatically by 
applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model, and then selecting the topology with 
superior log likelihood value. Sequences used for analysis are provided in Table S2 in File SI. 
doi:1 0.1 371/journal.pone.0089788.g001 
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Figure 2. Structural movements of RIG1 and MDA5 receptors during MD simulation. Porcupine plots showing structural movements of (Aj 
RIGl and (B) MDA5 receptors during IVID simulations. The structural movements are shown for both RNA-bound and free proteins. The docked dsRNA 
molecule has been colored in violet color. Arrow heads indicate direction of motion and length of the arrows specifies the extent of displacement. 
The figure also shows the schematic representation of intermolecular interactions (C) RIGl and (D) MDA5 with dsRNA after MD simulations. 
doi:1 0.1 371 /journal.pone.0089788.g002 



The intermolecular contacts of bufTalo RIG 1 and MDA5 with 
dsRNA was identified and compared with those of human 
(Figure 2C-D; Table 3). We found that buffalo and human share 
a similar pattern of dsRNA recognition mechanism. In buffalo and 
human RIG 1 , positively charged residues were primarily involved 
in dsRNA recognition. Additionally, hydrophobic and polar 
residues were found to interact with dsRNA. In contrast, polar 
and hydrophobic residues of MDA5 were mostly involved in 
dsRNA binding before MD simulation. However, the number of 
Arginine around dsRNA increased considerably after MD 
simulation (Table 3). On the contrary, MDA5 displayed greater 
number of Lysine around dsRNA, suggesting positively charged 
residues are indispensible for viral RNA recognition by buffalo 
RIGl and MDA5. 

The overall binding free energy (AGbind) of MDA5-dsRNA 
complex was found to be greater than that of the RIG 1 — dsRNA 
complex (Table 4). The affinity of dsRNA towards the receptors 
differs only at the electrostatic (AG,.„ui) and polar contributions 
(AGpoiar). The AG„,ui and AGp„i,,^ of MDA5-dsRNA complex 
were calculated to be greater than those of RIG 1 — dsRNA 
complex. Hence, it appears that polar and electrostatic interac- 
tions play dominant role for dsRNA recognition by RIGl and 
MDA5. 

Constitutive Expression of RLR Genes in Different Tissues 

Basal expression of the three buffalo RLR genes (RIGl, MDA5, 
and LGP-2) was studied in liver, spleen, intestine, kidney, and 
PBMC using qRT-PCR with RPS-18 (40S ribosomal protein SI 8) 
as the internal reference gene. It was found that the RLR genes 
were ubiquitously expressed in every tissue tested (Figure 3). The 



expression levels of RLR genes were comparatively higher in 
kidney and intestine than in spleen and PBMC. The ubiquitous 
nature of RLR expression of was also reflected at protein level. 
Immunohistochemical localization of RLR proteins were detected 
in primary immune organs (spleen and tonsil) as well as in thyroid, 
kidney, which are not directiy involved in primary immune 
response (Figure 4). Marked staining can be seen for RIGl and 
MDA5 across the red pulp area of spleen, follicular epithelial cells 
and parafollicular cells of thyroid tissue, renal tubular epithelial 
cells, mucosal and sub-mucosal layers of intestine. A lower 
expression of RIG 1 and MD A5 was observed in crypts of tonsil 
and in hepatocytes. Negative control sections with an isotypic IgG 
antibody were run in parallel; no staining was observed in any of 
these control sections. The primary antibodies used for RIGl and 
MDA5 were specific, as revealed from western blot analysis 
(Figure S5 in File SI). 

Responses of RLR and Associated Genes to poly l:C 

Following poly I:C exposure, mRNA expression pattern of RLR 
genes, downstream adapter genes (MAVS, IRF-3 and IRF-7) and 
the effector genes (IFN-B and ISGs) were analyzed in foetal 
fibroblast cells. Fibroblast cells provide an obvious advantage for 
studying RLR mediated antiviral response, because TLR .system 
for viral recognition is dispensable in this cell type [30,31] Primary 
culture was chosen over established immortahzed cell lines because 
the later may have multiple mutations that can affect the 
permissiveness and responsiveness to viral infections [32]. All 
experiments were carried out on fibroblasts with less than 10 
passages. 
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Figure 3. Constitutive mRNA expression levels of buffalo RLR {R/G7, MDAS and LGP2) genes in different tissues. Relative abundance of 
mRNA was measured by qRT-PCR using RPS-18 as housekeeping control. 
doi:1 0.1 371 /journal. pone.0089788.g003 



Expression Patterns of RLR Genes 

As shown in Figure 5, RLR genes showed no change in 
expression after 1 h post treatment, but RIGl and MDA5 genes 
were up-regulated (p<0.01) after 4 h post treatment. Maximum 
expression of tliese genes was seen during 6 h post treatment. 
Afterwards, the expression oi RIGl was sustained at a higher level 
until 48 h. The level of MDA5 dropped transiently during 1 2 h. 



but increased again during subsequent intervals of experiment. 
The expression pattern of LGP2 followed a slightly different 
pattern (Figure 5C). The mRNA expression level of LGP2 was 
significantly higher (p<0.01) during 6 h and was maximum at 
24 h after poly I:C treatment. It appeared that the amplitude of 
response was higher for MDA5, while RIGl maintained a sustained 
response following poly I:C treatment. 



Table 3. Comparative intermolecular interaction of RIGl and MDA5 (hunnan and buffalo) with dsRNA. 
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Table 4. MM/PBSA binding free energies (kJ/mol) of RIGl/ 
dsRNA and MDA5/dsRNA complexes. 
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AGcoui = Electrostatic energy. 

AGps = Polar solvation energy. 
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AGr,ps = Nonpolar solvation energy. 

AG„onpoiar = Nonpolar contribution (AG„ps+AG„dw). 

AGbind- Overall binding energy. 

doi:10.1371/journal.pone.0089788.t004 

Expression Patterns of MAVS, IRF-3, and IRF-7 

The expression pattern of MAVS showed an unusual trend 
following poly I:C treatment (Figure 6A). The mRNA expression 
level was significantly lower (p<0.05) during 4 h post-treatment. 
The expression level increased during subsequent hours and was 
significandy higher (p<0.01) during 12-24 h post-treatment. 
Again during 48 h the expression dropped to basal level. The 
induction of both IRF-3 and IRF-7 following poly I:C treatment 
was evident in the present experiment (Figure 6B-C). The mRNA 
levels of both IRF-3 and IRF-7 were significantly higher during 
6 h post treatment, though amplitude of response was much 
higher for IRF-7. The expression level showed a marked increase 
(p<0.01) during 24 h for both IRF-3 and IRF-7. In general, the 
expression pattern of IRF-3 and IRF-7 followed a similar trend 
until 24 h, after which the level of IRF-3 was constant until 48 h, 
but the level of IRF-7 dropped close to basal level. 



Expression Patterns of IFN-(3 and ISGs 

Poly I:C was found to induce IFN-P response in bovine foetal 
fibroblast cells. As shown in Figure 7A, IFN-P expression started 
increasing after 4 h post-treatment. Expression of IFN-P mRNA 
was highest during 6 h, higher level was maintained up to 12 h, 
and finally returned to basal level at 24 h post treatment. 
Induction of Type-I IFN response leads to transcriptional 
induction of >300 IFN-stimulated genes (ISGs) that have been 
shown to function as antiviral effectors. This was evident in the 
present experiment as well. As early as 4 h post-treatment, the 
transcript level of all ISGs (ISG-15, ISG-54 and ISG-56) increased 
(Figure 7B-D). At 6 h post treatment, a sharp increase in the 
expression level of all three ISGs was observed, the level of ISG-15 
being the maximum. It was found that during subsequent time 
period the levels of all ISGs were sustained at higher concentra- 
tion. 

Discussion 

Over the last decade, studies have highlighted the biological role 
and importance of RLRs in antiviral immunity [7]. However, 
compared to other mammalian, avian, amphibian, and fish [33— 
35] species, bovine RLRs have been poorly studied. Here, we 
reported the evidence of functional RIGl and MDA5 genes in 
buffalo genome. Domain analysis revealed that buffalo RIGl and 
MDA5 receptors have common domain architecture as in other 
mammals. It is noteworthy that only MDA5, but no RIGl 
ortholog, is present in all teleosts, the first chordates to be evolved. 
This indicates that MDA5 might precede RIG 1 in the evolution- 
ary history. It has been postulated that the functional domains 
present in MDA5 and RIGl have evolved independently, and 
most probably, not by a simple gene-duplication event [36]. The 
helicase domains might have been evolved from a eukaryotic 
ancestor gene, common to DICER and DEAD/DEAH box RNA 
helicase-I of MDA5 and RIGl, but the emergence of distinct 
CARD domains might have happened by domain grafting 
[36,37]. However, identification of RIGl-like genes in the 
genomes of sea anemone [Nematostella vectensis) and sea urchin 
[38] suggests that RIGl might have emerged from invertebrates, 
but gradually got lost in some fish lineages [37] . While the origin of 
RLR genes remains inconclusive, our sequence analysis indicated 




Figure 4. Constitutive protein expression levels of buffalo RIGl and iVIDAS in different tissues. Imuunohistochemical localization of RIGl 
and IVIDAS were detected in tfiin tissue sections using specific antibodies witfi DAB (3,3'-Diaminobenzidine) as substrate. Counterstaining was 
performed with Mayer's haematoxylin. 
doi:1 0.1 371/journal.pone.0089788.g004 
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Figure 5. Poly l:C induced expression levels of RLR genes in buffalo foetal fibroblast cells over different time intervals. Relative fold 
change in mRNA of (A) R1G1 (B) iVlDAS and (C) LGP2 over control has been shown. 
doi:1 0.1 371 /journal.pone.0089788.g005 



that the DEAD/DEAH box RNA hehcase (HEL-1 and HEL-2) of 
these two genes may share a common ancestor as revealed by the 
conservation of critical binding site residues viz. ATP binding site 
[G-(C/S)-G-K-T and D-Q-R-R], Mg++ binding site [D-E-C-H], 
and nucleotide binding site [K-T-R-(A/Q)-(P/H)-(I/V/S/N)-T- 
(S/T)-V] . Despite high sequence identity at RNA binding sites (E- 
Q/H-I/V-K) of RIG 1 and MDA5, substantial differences were 
observed within HEL-2i of both the proteins. However, the 
residues of HEL-2i domain were conserved in both RIGl and 
MDA5 across other mammalian species considered for this study. 
Our result suggests that irrespective of their point of phylogenetic 
diversion, critical binding site residues of RIG 1 and MDA5 are 
extremely conserved. This indicates that RNA binding pattern, a 
critical step for inducing antiviral innate immune responses, must 
be similar throughout the mammalian species. 

High resolution crystal structures of RIG 1 and MDA5 with or 
without viral RNA have provided a better understanding of RNA- 
induced receptor activation [5,10,39-42]. However, these studies 
were carried out on human, mouse, or duck models, but 
information on RNA mediated receptor activation in bovine is 
sparse. Our simulation studies elucidated a network of interactions 
between dsRNA and four domains of buffalo RIGl and MDA5. 
This mode of interaction was consistent with those observed in the 
crystallographic structures [5-6,39]. In both the receptors, 
positively charged residues were involved in providing affinity to 
dsRNA, as observed by the binding free energy calculation. 
Analysis of binding free energy terms indicated that MDA5 binds 
dsRNA with a higher affinity than RIGl does. We found the 
number of positively charged amino acids that interacted with 
dsRNA were higher in MDA5 than RIGl. This could be the 
reason of difference in net binding free energy and affinity of 
dsRNA-receptor complexes. 



The constitutive expression levels of the RLR genes were not 
restricted to the immune organs. Instead, organs like kidney and 
intestine showed a higher level of constitutive expression. This 
ubiquitous expression pattern of RLRs has also been observed in 
human, mice, and fishes [43,44]. Higher expression level in 
intestine is consistent with these receptors' role in maintaining the 
intestinal epithelial barrier that prevents entry of many viruses. 
Taken together, these findings suggest that many cell-types have 
the potential for RLR mediated response against viral invasion. 

The inductive expression patterns of RIGl and MDA5 in 
response to poly I:C were different. The induction of MDA5 
started early, amplitude of expression was higher, but expression 
level varied over time. However, RIGl was induced later, but its 
expression level was sustained during different hours of experi- 
ment. The higher amplitude of MDA5 induction could be 
attributed to its greater affinity for dsRNA than RIG 1 , as revealed 
from our MD simulation studies. It seems that LGP2 negatively 
regulated both RIGl and MDA5, as its expression pattern 
followed a reverse trend compared to RIGl and MDA5. The role 
of LGP2 in RLR mediated antiviral responses remains contro- 
versial as it has been reported both as a suppressor [19,20] as well 
as a positive inducer [21,22] of RIGl and MDA5 elicited 
signaling. LGP2 lacks N-terminal CARD regions, essential for 
MAVS mediated signaling, and keeps a check on RIGl and 
MDA5 mediated overproduction of IFNs [45]. Our results 
suggested that while RIGl, MDA5, and LGP2 are simultaneously 
involved in viral RNA recognition and the role of LGP2 might be 
regulatory in nature. 

After sensing dsRNA, RIGl and MDA5 interact with 
CARD domain of MAVS, leading to downstream activation of 
IRF-3, IRF-7 and IFN-(5 production [2,46]. However, we found 
the mRNA expression of MAVS was significantly (p<.05) 
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Figure 6. Poly l:C induced expression levels of RLR adapter genes in buffalo foetal fibroblast cells over different time intervals. 

Relative fold change in mRNA of (A) MAVS (B) IRF-3 and (C) IRF-7 over control has been shown. 
doi:1 0.1 371/journal.pone.0089788.g006 
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Figure 7. Poly l:C induced expression levels of IFN-p and ISGs in buffalo foetal fibroblast cells over different time intervals. Relative 
fold change in mRNA of (A) IFN-|3 (B) ISG15 (C) ISG54 and (D) ISG56 over control has been shown. 
doi:1 0.1 371 /journal.pone.0089788.g007 



downregulated during 4 h. RLR activation can induce a selective 
proteasomal degradation of MAVS, whicli is essential for 
downstream signaling [47]. Also, poly I:C can significandy 
downregulate MAVS expression in early hours of induction by 
rapidly decaying the MAVS mRNA [48]. Therefore, the initial 
downregulation of MAVS could be poly I:C mediated degradation 
of mRNA. Subsequent increase may be a consequence of cells' 
demand of MAVS for RLR signaling. MAVS activation causes 
homo and hetero dimerization of IRF3 and IRF7, which induce 
type I IFNs. IRF3 is constitutively expressed and involved in 
immediate antiviral responses, while IRF7 is expressed in a limited 
fashion and is induced after IFNP production is initiated [49]. 
Consistent with these findings, the present study also showed 
delayed induction of IRFs. Expression of IRF7 was temporal with 
high amplitude, but IRF-3 consistently showed a higher expression 
level after induction. 

Following poly I:C treatment, buffalo foetal fibroblast cells 
triggered an elevated IFN-|3 expression and transcription of IFN- 
stimulated genes (ISGs). The expression pattern of IFN-P was 
consistent with previous studies [50]. It is worth noting that despite 
negligible expression level of IFN-P during 24—48 h, the expres- 
sion levels of ISG15 and ISG54 continued to be high. This could 
be explained as activated homo- and heterodimers of IRF3 and 
IRF7 are able to bind IFN-stimulated response elements (ISRE) at 



the promoter regions of ISGs and can direct IFN-independent 
activation of ISGs [51]. Among ISGs, ISG56 can negatively 
regulate ISG15, ISG56, IRF3 and IFN-|3 by interacting with die 
adapter protein MITA to disrupt virus-induced IRF3 activation 
and cellular antiviral responses [52] . In agreement to this finding, 
the expression pattern of ISGs in our study suggests that ISG56 
could be a negative regulator of ISG15 and ISG54, besides 
inhibiting the expression of IFN-fi. 

Conclusion 

In conclusion, functional RIGl and MDA5 genes of buffalo 
were found to foUow similar pattern of domain organization as in 
other mammals. MD simulation studies suggested that despite 
striking homology at binding site residues, RIGl and MDA5 
receptors have different affinity towards dsRNA. Buffalo RLR 
genes are ubiquitously expressed in different tissues, irrespective of 
function of the tissue as immune organ. Poly I:C can induce an 
antiviral response in buffalo foetal fibroblast cells through RLR 
mediated pathway. 

Supporting information 

File SI Supporting Text, Tables, and Figures. Text SI: 

Methodology for isolation and culture of buffalo foetal fibroblast 
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cells. Text S2: Procxdures for model construction and simulation 
analysis. Text S3: Deduced amino acid sequence of buHalo RIGl 
and MDA5. Table SI: PGR cycling parameters used for the 
amplification of RLR genes. Table S2: Sequences (NCBI 
Accession numbers) used for phylogenetic analysis. Table S3: 
Comparative analysis of the stereo-chemical parameters of RIGl, 
MDA5 and their corresponding templates. Figure SI: Agarose 
gel electrophoresis of amplified products of different fragments of 
buffalo RIGl and MDA5 genes. Figure S2: Overall tertiary 
structure of helicase and C-terminal domains of (A) RIGl and (B) 
MDA5 proteins of buffalo. The helicases and C terminal domains 
of buffalo RIG 1 and MDA.5 were modeled based on the X-Ray 
crystallographic structures of human RIG 1 and MDA5 (termed as 
templates). The pair- wise alignment of buffalo and human 
sequences showed 82 and 84 percent sequence identities for 
RIGl and MDA5, respectively. Due to this striking homology 
between target and template proteins, the built models retained all 
the key structural features of the templates that included HEL- 1 , 
HEL-2, HEL-2i, and the C-terminal domains. Figure S3: 
Ramachandran plot of buffalo (A) RIGl and (B) MDA5 receptors. 
Stereochemical qualities of modeled proteins were found to be 
highly comparable to those of templates, indicating quality of the 
models were reasonably good to carr^- out further studies. Figure 
S4: Stability parameters of RIGl and MDA5 receptors as a 
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